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Abstract 

The first outbreak of dengue occurred in Madeira Island on 2012, featuring one virus 
serotype. Aedes aegypti was the vector of the disease and it is unlikely that it will be elimi¬ 
nated from the island. Therefore, a new outbreak of dengue fever can occur and, if it happens, 
risk to the population increases if two serotypes coexist. In this paper, mathematical model¬ 
ing and numerical simulations are carried out to forecast what may happen in Madeira Island 
in such scenario. 
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1 Introduction 

Dengue is the most rapidly spreading mosquito-borne viral disease in the world |46l I48j . In the 
last decades, dengue’s incidence has grown 30-fold and is now endemic throughout tropical and 
subtropical regions, present in more than 100 countries. According to World Health Organization 
(WHO), over 40% of world’s population is at risk 0131]. In Europe, in the last decade, some 
outbreaks have occurred [H] [23l [36] . 

Dengue is transmitted from an infected human to an Aedes mosquito, mainly Aedes aegypti^ 
and then from the mosquito on to other humans. The virus, which is circulating in the blood of 
the infected subject, is ingested by a female mosquito when it bites. After the incubation period, 
the infectious mosquito can transmit the virus on to other humans in subsequent bites |47) . The 
mosquito life cycle has four distinct stages: egg, larva, pupa and adult. The first three stages take 
place in the water, while the adult stage occurs in the air. 

*This is a preprint of a paper whose final and definite form is in Operations Research For Health Care, ISSN 2211- 
6923. Paper submitted 21/Nov/2014; revised 17/May/2015 and 07/Jul/2015; accepted for publication 14/Jul/2015. 
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Dengue can be caused by four different serotypes, named DEN-1,2,3,4. It is well documented 
in the literature that there are differences among the serotypes (see, e.g., [3]), but it is still to 
be determined which clinical characteristics correspond to which serotype [3]. However, several 
reports have indicated that serotypes DEN-2 and DEN-3 may cause more severe disease than 
the serotypes DEN-1 and DEN-4 (see, e.g., [331 Infection by one serotype provides lifelong 
immunity against that virus and temporary immunity against the other three. After a cross¬ 
immunity period of 2 to 3 months, a subject can be infected with a different serotype [nms]. 
The existence of four different viruses causes a wide clinical spectrum that includes asymptomatic 
cases, classic cases of dengue fever and more severe cases known as Dengue Hemorrhagic Fever 
(DHF). There is good evidence that a subsequent infection by other serotype increases the risk of 
developing severe dengue BEE]. 

As there is no vaccine or specific treatment for dengue, prevention is very important. To prevent 
the dengue virus transmission, the mosquito vector must be controlled to avoid reproduction 
and consequent bites. Control of vector is mainly achieved by eliminating the places where the 
mosquito lays their eggs, such as artificial water containers. As adult mosquitoes bite inside or 
outside houses, during the day or at night when the lights are on, it is also essential for humans 
to apply repellent on skin and also to use personal household protection like window screens, 
insecticide treated materials and long-sleeved clothes BEidn]. 

Rich literature investigating the dengue disease from the mathematical point of view exist. 
Some run simulations on specific regions of the globe HOIIMIISIIISIISS] and others make studies 
using the effects of the disease spread considering biological issues m Hi 1301131 [35]. In this 
paper we aim to study a preview scenario for Madeira Island in the case of a second outbreak 
with two serotypes circulating. This is important because DHF is 15 to 80 times more likely in 
secondary infections, being positively associated with pre-existing dengue virus antibodies m- 

In the Portuguese island of Madeira, the first outbreak of dengue fever occurred between 
October 2012 and February 2013. This outbreak was caused by only one serotype, DEN-1 virus, 
and it was declared controlled March 12, 2013. The vector was the mosquito Aedes aegypti, 
detected in Madeira for the first time in 2005. The mosquito is still not eradicated from Madeira 
and the Portuguese Ministry for Health consider that it can hardly be miniin]. The island has 
an intensive trade with areas where the disease is already endemic and the number of tourists has 
been increasing. These two factors, beneficial to the local economy, can also lead to the appearance 
of a new serotype in the island and an outbreak can occur at any moment |38j . 

The main motivation of this work is to analyze what may happen in Madeira Island in a 
dengue outbreak if two serotypes of dengue coexist. Mathematical modeling is used to simulate 
this hypothetical outbreak, caused by two different serotypes, considering the use of insecticide as 
a control measure. The article is organized as follows: Sectionprovides a brief description of the 
first outbreak in Madeira and in Section [3] we formulate the mathematical model to describe the 
interaction between human and vectors. Numerical experiments are presented in Section and 
conclusions ensue in Section |5l 


2 Madeira’s 2012 outbreak 

The outbreak affected mainly Funchal, the main city of the island, and lasted for 21 weeks, ending 
late February 2013. By that time, 2171 probable cases of dengue fever had been reported and, of 
these, 1084 were confirmed in the laboratory. In Figurej^the reported cases per week are presented. 
No deaths or cases of DHF were reported gOlllS]. In order to control the mosquito vector, local 
authorities took a whole range of measures, including the usage of a massive quantity of insecticide, 
starting even before the outbreak. However, Aedes aegypti has shown a high level of resistance 
against the insecticide [38]. As a response, educational campaigns were conducted. Media coverage 
in newspapers, television, radio, the internet and through flyers was a powerful tool. In particular, 
the population was informed about the importance of using repellent and covering arms and legs. 
Artificial water containers, such as cans, jars, barrels, plastic containers and broken bottles were 
requested to be removed. In addition, a medical appointment, exclusively dedicated to dengue 
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Figure 1: Notified dengue fever cases in Madeira per week, from October 2012 to February 2013 
(Source: Institute de Administragao da Saiide e Assuntos Sociais, Regiao Autonoma da Madeira) 


fever, was implemented in Bom Jesus Health Centre. Nowadays, local authorities continue to 
monitor Aedes aegypti entomological activity through a program of traps monitoring |121140] . 


3 The mathematical model 


In order to describe what can happen if two virus serotypes coexist, a mathematical model based 
on [36] is here proposed. The novelty is the presence of a second serotype of dengue and the 
consequent presence of a new set of variables: 

Sh — susceptible human (who can contract the disease); 

Iih — human infected by serotype 1 (who can transmit the serotype 1); 

Rih — human resistant to serotype I (infected by serotype 1 and recovered); 

Ijh — human infected by serotype j (who can transmit the serotype j); 

Rjh — human resistant to serotype j (infected by serotype j and recovered); 

Iijh — human infected first by serotype 1 and after by serotype j; 

Rijh — human first resistant to serotype 1 and after to serotype j; 

Ijih — human first infected by serotype j and after by serotype 1; 

Rjih — human first resistant to serotype j and after to serotype 1; 

Sm — susceptible mosquitoes; 

Am — aquatic phase (stages that take place in water: egg, larva and pupa); 
hm — mosquitoes infected by serotype 1 (that can transmit the disease); 

Ijm — mosquitoes infected by serotype j (that can transmit the disease). 


The index 1 refers to the virus serotype DEN-1, the responsible for the first outbreak in Madeira 
Island. The new index j, j G {2, 3,4}, refers to different virus serotypes (DEN-2, 3 or 4). 

The first nine state-variables are related to human population while the last four are concerned 
with the female mosquito. These compartments are mutually-exclusive m- It should be noted 
that there is no variable to represent resistant mosquitoes or to represent mosquitoes infected by 
a second serotype. Indeed, given that mosquitoes have a very short lifespan, there is no time for 
them to recover from an infection. Thus, once infected, the mosquito remains infected until its 

death (SET]. 

The application of insecticides is the most common control measure. In order to simplify the 
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model, insecticide was the only control measure considered. More precisely, we consider Cm(t) to 
be the proportion of insecticide applied at time t, with 0 < Cm(t) < 1- 
The system of differential equations for the human population is: 


^ = fihNh - + /i,,) Sh 

at \ Nh Nh J 

= BPimh^^Sh - {rjih + ^^h) hh 
at Nh 

= BlSjmh^Sh - (jljh + ^lh) Ijh 

at i\h 
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dRjh j ( Q ^Im \ 
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The differential equations are subject to the initial conditions 


Sh{t)) — SfiQ^ — diHO: djh{t}') — djfiO^ dijfiiO^ — djifi{0') — djifiO^ 

— BlhO: — BjhO: — BljhO^ Bjlh{t}') — BjlhO: 

^m(O) — -^mOi — SmOi .^lm(O) — — .fjmO- 

A wide range of parameters is used to build the model, e.g., by Nh we denote the total human 
population in Madeira; mortality is equal to the inverse of lifespan with a constant rate to 
human population and /i^ for the mosquito population; an infected mosquito can transmit dengue 
to a susceptible human with a probability Pmh per bite while an infected human transmits the 
disease to a susceptible mosquito with a probability jdhm per bite. Table presents the model 
parameters, their meaning and their values. An explanation is, however, necessary. 

There is evidence that a secondary infection increases the risk of developing DHF, which can 
be explained by the ADE phenomenon {Antibody-Dependent Enhancement) |44| . Recovering from 
one serotype of dengue, the immune system responds by producing antibodies to the virus, gaining 
lifelong immunity against that serotype. However, when a new serotype of dengue is contracted, 
the same antibodies that protect against the previous serotype, facilitate the entry of the new 
virus into host cells, enhancing the infection, with an increase in viral production. Viral loads are 
associated with transmissibility, and several researchers believe that subjects are more infectious 
in the second infection than during the hrst, and consequently the transmission rate increases [5]. 
In this paper we use p,dhf to represent death from DHF and cr to represent the ADE phenomenon 
[mni. Note that the total population is not constant in our model because we are considering a 
probability of death from DHF different from zero. This is in contrast with [36]. Our approach to 
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the multistrain model consists of using different parameter values for different strains: rjih ^ rjjh 
and/or Pimh ^ Pjmh- For simplicity of notation, in Table we omit the specific strain 1 or j, 
JG {2,3,4}. 


Table 1: Parameters of the epidemiological model Q-© 


Para¬ 

meter 

Description 

Range of values 
in the literature 

Used 

values 

Source 

Nh 

total population 


112000 

m 

B 

average daily biting (per day) 


1/3 

m 

i^mh 

transmission probability from Im (per bite) 

[0.25, 0.33] 

0.25 

m 

j^hm 

transmission probability from Ih (per bite) 

[0.25, 0.33] 

0.25 

m 


average lifespan of humans (in days) 


79 X 365 

m 

i/% 

average viremic period (in days) 

[4, 15] 

7 

M 

i-^m 

average lifespan of adult mosquitoes (in days) 

[8, 45] 

15 

iminiiisi 


number of eggs at each deposit per capita (per day) 


6 

m 


natural mortality of larvae (per day) 


0.2363 

a 

VA 

maturation rate from larvae to adult (per day) 

[1/11, 1/7] 

1/9 

m 

k 

number of larvae per human 


0.9 

[niiig 

a 

ADE phenomenon 

[0, 5] 

0,5; 1.1; 2.5 

m 

k-dhf 

probability of death from DHF 

[0.01, 0.1] 

0.02 

lillTj 


In our model we assume that: 

• Infected mosquitoes do not transmit the virus onto their eggs; 

• The population is homogeneous; 

• There is no immigration or emigration; 

• There is no seasonality; 

• There is an homogeneous mixing of human and mosquito population, having the vector an 
equal probability to bite any human; 

• Humans and mosquitoes are assumed to be born susceptible. 

The scheme of the model 0-(Il is shown in Figure An arrow pointing into a compartment 
is associated with a positive member of the corresponding differential equation while an arrow 
pointing out of the compartment represents a negative member of the equation. 

Comparing our model with the previous one [36] . this one is more complex. The coexistence 
of two serotypes of dengue leads to 13 state variables, more than doubling the 6 state variables in 
the previous model |36| . As a consequence, the interaction between human and vector has become 
more complex and the number of arrows in Figure]^ has increased with respect to the scheme of 
[36]. It is also noted that there are no equilibrium points and, therefore, the basic reproduction 
number based on |42| cannot be computed. This is also in contrast with |36| . 

The impact of ADE on the model can be profound: single strain models usually cannot have 
isolated equilibria and show fluctuations without external seasonal drives or noise, while multi¬ 
strain models with ADE phenomenon may present oscillations |5]. 


4 Computational experiments 

The software used in our simulations was Matlab with the routine ode45. This solver is based on 
an explicit (4,5) Runge-Kutta formula, the Dormand-Prince pair. The numerical solver ode45 
combines fourth and fifth order methods, both of which are similar to the classical fourth-order 
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Figure 2: The scheme of the dengue model Q-© with two serotypes 


Runge-Kutta method. These vary the step size, choosing it at each step in an attempt to achieve 
the desired accuracy. 

Our computational experiments simulate a hypothetical outbreak in Madeira Island caused by 
two different serotypes 1 and j, j S {2,3,4}. We considered that the initial values related with 
serotype DEN-1 are the final values of Madeira’s outbreak of 2012. Funchal has a total population 
of 112 000. Of these, 2171 were infected. Serotype DEN-j is new, so = 0. As the other 

initial values are unknown, we presumed the following values: 

• lih{0) = C and ljh{0) = C; ? G {10,50,100} (hypothetical values, assuming that with these 
infected numbers, health authorities could act and take measures to control the disease); 

• lijh{0) = C) C € {10,50,100} (individuals who have recovered from the outbreak and are 
now infected by serotype j). 

As the first outbreak caused no deaths, all the infected have recovered from serotype 1. However, 
^ of these are now infected by serotype j, so 

Rm(0) = 2171-^, i?,;,(0) = 0, i?i,^(0) = 0, i?,ift(0) = 0. 

Thus, the initial value for susceptible humans is Sh{0) = 109829 — 2^. We considered different 
virulence and transmission intensities. Two scenarios are possible: or j € {2,3}, and i/r]jh and 
(3jmh are not lower than those of DEN-1 (more aggressive case); or j = 4, and l/rjjh and Pjmh 
do not have higher values than those of DEN-1 (less aggressive case). In order to evaluate the 
influence of the (3 and <j parameters and the insecticide control Cm, some numerical tests were 
carried out. At a first step, we fixed /3 = 0.25 and a = 1.1 and tested the model varying the 
control Cm- The corresponding results are shown in Figurej^ It should be noted that in Figuresj^ 
EE and 1^ the scale of the infected human axis is smaller in secondary infections than in hrst 
infections. This was essential due to the following reasons: 

• in the first infection people are naive to the disease; 
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(a) Iih (infected by DEN-1) (b) Ijh (infected by DEN-4) 




(c) Iijh (infected by DEN-1, then by DEN-4) (d) Ijih (infected by DEN-4, then by DEN-1) 

Figure 3: Infected individuals varying the control Cm in case of two dengue serotypes: DEN-1 and 
DEN-4 (results for 1 year with = 0.02, a = 1.1, = Pihm = 0.25, = Pjhm = 0.25, 

Vih = 1/7, rjjh = 1/5, and hho = Ijho = Iijho = 10) 


• there are a substantial number of resistant individuals; 

• there is a partial and transient immunity from the first infection to the second one, which 
delays second infections. 

Figure|^shows the fluctuation of the number of infected human with the application of different 
proportions of insecticide. It is patent that insecticide is a very efficient measure. Even with small 
quantities, its influence is huge. Note that with Cm = 0.02, in the considered time interval, the 
number of infected is tending to zero. 

It is clear that should a new outbreak of dengue with two different serotypes occur, health 
authorities should pay particular attention to the more aggressive serotype: if the new serotype 
is DEN-2/3, then the priority is to take care of people infected by DEN-2/3 (see Figure]^; if the 
new serotype is DEN-4, then health authorities should pay particular attention to DEN-1 (see 
Figure]^. In all cases, the control has an important role on the number of infected people and 
“doing nothing” should not be an option to health authorities. 

One important information is the evolution of the human population Nh- Figure [4a| presents 
this evolution considering (3 = 0.25 and a — 1.1 for different controls. Using no insecticide, 450 
deaths are expected at the end of 154 days. However, the application of a small quantity of 
insecticide reduces this number to 50. 

We now increase j3 to 0.33 keeping cr = 1.1 and varying the control Cm to study the influence 
of the transmission probability (see Figure [^. As expected, the number of infected people rises. 
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(a) DEN-1 and DEN-4 
(^l^jmh ~ (^jhm ~ 0.25, ~ 1/^) 


(b) DEN-1 and DEN-j, j € {2,3} 
{Pjmh = Pjhm = 0.33, r)jh = 1/9) 


Figure 4: Total population Nh varying the control Cm (results for 1 year with fidhf = 0.02, cr = 1.1, 
Plmh = Plhm = 0.25, r]ih = 1/7, /l/t(0) = IjhiO) = hjhiO) = 10) 




(a) Iih (infected by DEN-1) (b) Ijh (infected by DEN-/) 




(c) I\jh (infected by DEN-1, then by DEN-/) (d) Ijih (infected by DEN-/, then by DEN-1) 

Figure 5: Infected individuals varying the control Cm in case of two serotypes: DEN-1 and DEN-/, 
/ e {2, 3} (results for 1 year with = 0.02, a = 1.1, fiimh = Pikm = 0.25, I3jmh = (3jhm = 0.33, 
Vlh i/7, ^jh I/O, and I\hQ ^jho ^IjhO 10) 
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Time (days) Time (days) 


(c) Iljh (d) Ijlh 

Figure 6: Infected individuals varying cr, in case of two serotypes DEN-1 and DEN-4 


having increased almost four-fold. It is thus concluded that /3 is a very sensitive parameter. As 
before, the impact of the control is crucial. 

Figure 4b reports the evolution of the total population Nh using = 0.33 and a = 1.1 
for different controls. In this figure the application of insecticide is even more evident: with no 
control, human population can suffer a significant decline. With = 0.02, the fluctuation of total 
population is almost non-existing. Comparing with Figure [4a| it is also noticeable the influence of 
p. Increasing the transmission probability from infected individuals, the population fall is bigger. 
Note that Nh starts decreasing at t = 40, about the same time the number of secondary infections 
become significant (Figures [^c) and (d)). 

In the second phase of the numerical experiments, we fixed /3 = 0.25 and = 0.01 and tested 
the model varying a (see Figure]^. In Figures]^ a) and[^b) the influence of the parameter a is 
residual. Although the ADE phenomenon is just related to a secondary infection, a slightly affects 
the number of infected mosquitos and, consequently, has an impact, even though small, on the 
first infection. For this reason, this small influence only appears a few days after the beginning of 
the outbreak. In Figures [^c) and|^d) the ADE phenomenon is considerable. Increasing tr from 
0.5 to 2.5 leads to a five-fold increase of second infections {hjh and Ijih)- 

Figures and show the results for P increased to 0.33 while keeping the control Cm = 0.01 
and varying a. This allows to learn about the influence of the transmission probability P (cr seems 
to be more sensitive with P at 0.33). On one hand, it is evident in Figure]^ the influence of ADE 
phenomenon. On the other hand, when compared to Figure]^ it is obvious the importance of the 
P parameters. 
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(a) DEN-1 and DEN-4 (b) DEN-1 and DEN-j, j G {2, 3} 

Figure 7: Total population Nh varying a, in case of two serotypes 




Time (days) Time (days) 


(Sr) Ilh (^) 




Time (days) Time (days) 


(c) Iljh 


(d) IjiH 


Figure 8: Infected individuals varying cr, in case of two serotypes DEN-1 and DEN-j, j G {2,3} 


10 

















































(a) DEN-1 and DEN-4 (b) DEN-1 and DEN-j, j G {2, 3} 

l/Sjmh = Pjhm = 0.25, rjjh = 1/5) Ifijmh = /3jhm = 0.33, rjjh = 1/9) 

Figure 9: Total population Nh{t) varying the initial conditions /i;i(0) = ljh{0) = hMo) = e, 
^ G {10,50,100}, in presence of two serotypes (results for 1 year with = 0.02, Cm = 0.01, 
o- = 1.1, Pimh = Pihm = 0.25, ijih = 1/7) 


Additionally, we also study how many initial infections of new serotypes would allow these 
serotypes to persist in the population if no control measures were to be taken. For this, we set the 
value Cm = 0 and considered different initial values of infected people by the new serotype. More 
precisely, at time zero we considered /i;i(0) = ljh{0) = hjhiO) = C with ^ G {10,50,100}. For 
both serotypes DEN-2/3 or DEN-4 the curves have almost the same shape, the small variations in 
value being the expected ones: the curves hhit), luhit), hh(t), hihit), Ruhit), Rihit), R4,ih{t), 
l 2 / 3 h{t), h/sihit), Rih{t), Ri 2 / 3 h{t), R 2 / 3 h{t), R 2 / 3 ih{t) increase continuously with the increase 
of the initial value the curve Nh{t) decrease continuously with the increase of the initial value 

see Figure]^ Curiously, as can be seen in Figure [TOl an exception in behavior occurs in the 
first month for the curve j G {2, 3}. 

We finalize by remarking that the scales of figures for scenario j G {2,3} (serotype more 
aggressive than DEN-1) and for scenario j = 4 (serotype less aggressive than DEN-1) are quite 
different. Figure [TT] illustrates well how much more aggressive DEN-2 or DEN-3 is than DEN-4. 

5 Conclusion 

According to the World Health Organization (WHO), at least 40% of the world population is 
at risk from dengue and several outbreaks have appeared recently in Europe. The first was 
recorded in Greece and consisted of two serotypes, Madeira’s outbreak was the second epidemics 
of dengue in Europe. Local authorities in Madeira continue to monitor Aedes aegypti entomological 
activity through a program of traps nmn]. Moreover, the report [T3] of the European Centre for 
Disease Prevention and Control (ECDC) concludes that “The important increase of co-circulation 
of different dengue virus serotypes worldwide is a potential source of its re-introduction to Madeira 
in the future” and “The possibility of sustainability of the virus during the winter season and 
subsequent reemergence when climatic conditions again favor Aedes aegypti mosquito activity 
exists”. Motivated by these, we proposed an epidemic model for dengue fever infection that 
considers the simultaneous spreading of two different serotypes throughout shared non-naive host 
population, with different virulence and transmission intensity. Its dynamics was discussed in 
the context of Madeira island, and the influence of the infected transmission probability (human 
and mosquito) analyzed. The Antibody-Dependent Enhancement (ADE), which results from a 
new infection with a particular serotype in an individual with acquired immunity to a different 
serotype, was also studied. The higher the value of ADE, the higher is the number of infected 
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Time (days) 

Figure 10: Infected individuals first by DEN-1 and after by DEN-j, j G {2,3} (results for 1 year 
with fj^dhf ~ 0.02, Cttt, = 0.01, (7 — 1.1, ^imh ~ Pihm ~ 0.25, i^jmh ~ ^jhm ~ 0.33, T/x/i = 1/^7 and 
Vjh = 1/9) 



Figure 11: Total population Nh {p.dhf = 0.02, Cm = 0.01, a = 1.1, jSimh = Pihm = 0.25, 
Vih = 1/7, and hho = Ijho = hjm = 10) for DEN-4 {(3jmh = Pjhm = 0.25, rj^h = 1/5) versus 
DEN-j, j e {2, 3} {/3jmh = Pjhm = 0.33, r]jh = 1/9) 


12 

















individuals. More important, while classical dengue fever causes negligible mortality, ADE leads 
to the development of a significantly more dangerous dengue. The important role of insecticide as 
a control measure, is also a relevant conclusion. 

This study is useful for strategies of health control on Madeira Island. For example, we 
used one control variable related with application of insecticide. It turns out that application 
of insecticide was done in Madeira intensively, mainly carbamate formulations (bendiocarb) and 
Bacillus thuringiensis israelensis (bti), which is a bacteria-derived toxin that is lethal when ingested 
by larvae Indoor spraying at the biting time of the Aedes aegypti mosquito (early morning 
or late afternoon) is also easily done by population (dwelling rooms were sprayed by inhabitants 
of Madeira Island at the end of the day with doors closed, e.g., bathrooms, kitchens and toilets). 
However, the model investigated here is sufficiently general to be applied in other contexts. 
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